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A  new  systematic  approach  to  state  of  charge  (SoC)  observer  design  for  battery  cells  is  presented.  It  is 
based  on  a  purely  data  driven  model  and  a  nonlinear  observer  constructed  from  it.  As  a  key  novelty,  a 
unified  and  generic  framework  for  model  generation  (i.e.  design  of  experiments,  nonlinear  model 
structure)  and  observer  parametrisation  is  presented.  An  integral  part  of  SoC  observers  in  hybrid  elec¬ 
trical  vehicles  is  a  dynamic  battery  model  which  describes  the  nonlinear  system  behaviour  of  the  cell 
terminal  voltage.  In  order  to  enable  the  application  of  the  proposed  concepts  for  any  type  of  battery 
chemistry,  a  data  based  modelling  approach  using  the  architecture  of  local  model  networks  (LMNs)  is 
proposed.  As  an  important  prerequisite,  optimal  model  based  experiment  design  ensures  proper  exci¬ 
tation  of  the  system  dynamics  while  the  desired  SoC  operating  range  is  covered.  For  SoC  estimation,  an 
augmented  state  space  representation  of  the  LMN  is  derived  and  the  nonlinear  observer  design  is  pre¬ 
sented.  In  particular,  the  use  of  a  fuzzy  observer  is  beneficial  in  combination  with  LMNs  since  the  local 
observers  are  time-invariant  which  greatly  reduces  the  complexity  of  the  global  estimator.  The  proposed 
concepts  are  validated  experimentally  by  means  of  a  Lithium  Ion  power  cell. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

This  work  describes  a  new  systematic  approach  to  experiment 
design,  nonlinear  system  identification  and  observer  design  for 
battery  cells  for  the  purpose  of  state  of  charge  (SoC)  estimation. 

The  correct  determination  of  the  state-of  charge  and  also  the 
related  state  of  power  is  an  important  task  in  hybrid  electrical  ve¬ 
hicles  (HEV).  In  this  context,  the  traction  battery  is  a  complex 
system  consisting  of  many  single  cells  and  an  especially  designed 
electrical  circuit.  Part  of  that  circuit  —  and  by  far  the  most  critical 
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due  to  its  complicated  calibration  —  is  the  battery  management 
system  (BMS).  This  supervising  control  unit  processes  the  actual 
status  of  the  cells  of  the  traction  battery  and  shares  the  computed 
dataset  with  the  rest  of  the  vehicle’s  controllers  (in  particular  with 
the  hybrid  control  unit,  HCU).  It  determines  the  major  operational 
tasks  of  the  traction  battery  as  are  available  power  and  energy,  high 
reliability  (long  cycle  and  long  calendar  life),  and  battery  safety 
under  any  circumstances  (i.e.  use  and  abuse  cases).  Achieving  these 
goals  is  widely  connected  to  the  quality  of  the  battery  management 
system  and  its  interaction  with  the  other  vehicle  control  units  [1,2], 
An  integral  part  of  the  BMS  is  a  mathematical  cell  model  which 
describes  the  nonlinear  dynamic  behaviour  of  the  cell  terminal 
voltage  in  dependance  of  the  charge/discharge  current.  Such  a 
model  allows  to  predict  the  nonlinear  system  dynamics  of  the 
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traction  battery  under  the  specific  loads  and  environmental  con¬ 
ditions  and  is  required  for  SoC  estimation  since  it  is  not  possible  to 
measure  the  SoC  directly.  Typical  modelling  approaches  in  the  field 
of  battery  cells  are: 

•  Electro-chemical  modelling:  Such  models  are  based  on  a 
detailed  electro-chemical  description  of  the  cell,  see  e.g.  Refs. 
[3—5],  The  major  disadvantages  of  this  approach  are  that  time 
efficient  parametrisation  as  well  as  real-time  application  is 
very  complex  or  even  not  possible,  c.f.  Refs.  [6,7]. 

•  Equivalent  circuit  models:  These  models  comprise  a  combina¬ 
tion  of  RC  circuits  in  series  with  an  internal  resistance  and  an 
ideal  voltage  source.  However,  a  single  equivalent  circuit  model 
cannot  describe  the  battery  operation  over  a  large  range  of  SoC 
and  temperature,  [7],  Consequently,  an  individual  para¬ 
metrisation  for  various  operating  conditions  (temperature, 
SoC)  is  required. 

•  Black/grey  box  techniques:  Data  based  modelling  approaches 
offer  a  versatile  structure  for  the  identification  of  nonlinear 
systems  while  the  real-time  application  is  easily  possible,  e.g. 
Ref.  [8], 

Recent  publications  have  addressed  modelling  structures  for 
HEV  calibration:  In  Ref.  [7]  a  phenomenological  model  based  on  an 
equivalent  circuit  model  with  varying  parameters  is  presented.  The 
model  parameters  are  optimised  based  on  measured  data  using  a 
genetic  algorithm.  Another  interesting  approach  using  state  space 
model  structures  with  additional  hysteresis  and  filter  states  is 
described  in  Ref.  [6]. 

A  widely  used  approach  for  SoC  estimation  of  batteries  is  the 
extended  Kalman  filter  (EKF)  in  combination  with  an  equivalent 
circuit  model,  see  e.g.  Refs.  [9—12].  The  online  estimation  of  the  SoC 
of  a  Lithium  Ion  cell  based  on  an  electro-chemical  model  can  be 
found  in  Ref.  [13],  A  simple  resistor-capacitor  battery  model  is  used 
in  Ref.  [14]  where  modelling  errors  caused  by  the  simple  model  are 
compensated  by  a  sliding  mode  observer.  In  Ref.  [8]  the  battery  is 
modelled  using  a  radial  basis  function  neural  network  and  the  state 
of  charge  estimation  based  on  the  EKF  is  presented.  In  order  to 
apply  the  EKF  to  the  proposed  nonlinear  model,  the  battery  model 
must  be  linearised  at  every  sampling  instance.  A  very  similar 
approach  can  be  found  in  Ref.  [15],  In  Ref.  [6]  the  modelling  and 
identification  of  a  state  space  structure  which  includes  terms  that 
describe  the  dynamic  contributions  due  to  open  circuit  voltage, 
ohmic  loss,  polarization  time  constants  is  described.  Again,  an  EKF 
is  used  to  estimate  the  SoC. 

In  this  paper,  a  systematic  and  integrated  approach  to  cell 
modelling  and  the  associated  nonlinear  observer  design  is  pre¬ 
sented  and  the  results  are  validated  by  means  of  real  measurement 
data  from  a  Lithium  Ion  cell.  However,  the  proposed  method  aims 
at  time  efficient  modelling  procedure  suitable  for  different  cell 
types/chemistries  and  the  subsequent  automated  observer  design 
based  on  the  nonlinear  cell  model.  Therefore,  a  data  based 
modelling  approach  is  chosen  due  to  the  applicability  for  various 
types  of  batteries.  Additionally,  the  proposed  nonlinear  model 
structure  allows  for  an  efficient  and  automated  nonlinear  observer 
design.  The  systematic  approach  to  SoC  observer  design  presented 
in  this  paper  can  be  subdivided  into  three  steps: 

•  Optimal  experiment  design:  Experiment  design  is  an  important 
prerequisite  for  data  based  modelling  and  has  the  goal  to  gain 
maximum  of  information  from  an  unknown  system  with 
minimum  experimentation  effort.  Especially  with  battery  cell 
modelling,  the  design  of  experiments  (DoE)  plays  an  important 
role  due  to  the  fact  that  the  SoC  excitation  (and  the  SoC  oper¬ 
ating  range)  directly  depends  on  the  excitation  signal  of  the  cell 


current.  In  this  work,  a  model  based  experiment  design 
methodology  for  battery  cells  is  presented,  where  the  genera¬ 
tion  of  the  excitation  signal  is  particularly  adjusted  to  the 
system  in  that  a  prior  process  model  is  available,  see  also  Ref. 
[16],  In  general,  the  target  of  model  based  DoE  is  that  param¬ 
eters  belonging  to  a  specific  model  structure  can  be  estimated 
from  measured  data  with  minimal  variance  [17].  In  addition, 
the  optimal  experiment  design  procedure  helps  to  reduce  time 
and  costs  for  cell  testing. 

•  Nonlinear  system  identification:  The  SoC  observer  design  is 
based  on  a  nonlinear  model  of  the  battery.  In  the  context  of 
data  based  modelling,  local  model  networks  (LMNs)  are  a 
widely  used  concept  in  automotive  industry,  see  e.g.  Refs.  [18— 
20].  These  models  interpolate  between  different  local  models, 
each  valid  in  a  certain  operating  regime,  see  e.g.  Refs.  [21,22]. 
They  offer  a  versatile  structure  for  the  identification  of 
nonlinear  dynamic  systems  and  the  incorporation  of  prior 
(physical)  knowledge  is  easily  possible  due  to  the  transparency 
of  the  LMN  architecture.  In  the  context  of  cell  modelling,  spe¬ 
cial  effects  such  as  relaxation  and  hysteresis  impose  an  addi¬ 
tional  challenge.  Thus,  the  based  on  Ref.  [23]  (and  Ref.  [24] 
respectively),  the  adaptation  and  extension  of  LMN  structures 
for  cell  identification  is  addressed  in  this  paper. 

•  Nonlinear  observer  design:  LMNs  offer  local  interpretability  such 
that  local  linear  battery  impedance  models  can  be  extracted  and 
thus  provide  a  basis  for  automated  nonlinear  observer  design.  In 
this  work,  fuzzy  observer  design  based  on  the  proposed  model 
structure  for  the  purpose  of  SoC  estimation  is  presented.  The 
fuzzy  observer  uses  standard  Kalman  filter  theory  for  each  local 
linear  state  space  model.  Linear  combinations  of  the  local  filters 
are  then  used  to  derive  a  global  filter  [25],  The  biggest  advantage 
compared  to  the  widely  used  extended  Kalman  filter  is  that  the 
local  observers  are  time-invariant  and  no  linearisation  at  every 
sampling  instance  is  required,  which  greatly  reduces  the 
computational  complexity  of  the  global  filter. 

The  remainder  of  this  paper  is  structured  as  follows:  In  Section 
2,  SoC  determination  in  general  and  the  SoC  observer  architecture 
is  addressed.  In  Section  3,  the  battery  cell  modelling  workflow 
including  model  based  experiment  design  is  presented.  Addition¬ 
ally,  the  architecture  of  the  LMN  and  the  training  algorithm  shortly 
reviewed.  Design  and  architecture  of  the  fuzzy  observer  are 
described  in  Section  4.  The  proposed  concepts  are  validated  by 
means  of  real  measurement  data  from  a  Lithium  Ion  cell. 

2.  Nonlinear  SoC  observer 

In  this  section,  SoC  determination  in  general  and  model  based 
SoC  estimation  using  different  observer  structures  are  discussed. 

Basically,  there  are  two  practical  approaches  how  to  determine 
the  SoC  in  a  mobile  battery  [26]:  One  is  the  relative  SoC  determi¬ 
nation  where  an  initial  SoC  is  adjusted  by  taking  the  time  integral  of 
the  current: 

SoC(t)  =  SoC0  +  j  ^Mlldr  (1) 

i  =  0 

where  SoC0  denotes  the  initial  SoC,  /(t)  the  instantaneous  cell  cur¬ 
rent,  Cn  the  nominal  cell  capacity  and  rji (/)  is  the  coulombic  effi¬ 
ciency.  The  drawbacks  of  this  method  are  that  the  initial  SoC  has  to 
be  known  and  that  the  relative  SoC  determination  becomes  unre¬ 
liable  when  operated  in  this  mode  for  an  extended  time.  Addi¬ 
tionally,  the  internal  resistance  of  the  battery  causes  thermal 
energy  losses  upon  (dis-)charging. 
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Fig.  1.  OCV  measurement  during  charging  and  discharging. 

Secondly,  it  is  possible  to  use  the  open  circuit  voltage  (OCV)  of  a 
battery  cell  for  an  absolute  SoC  estimation.  Since  the  open  circuit 
voltage  is  the  only  measurable  battery-intrinsic  variable,  this 
approach  is  more  precise  in  most  cases.  However,  it  is  not  possible 
to  directly  determine  the  SoC  using  the  OCV  when  the  unit  is  under 
operation,  due  to  a  hysteresis  between  the  charging  and  dis¬ 
charging  voltage  profile  which  depends  on  the  load  current.  In  Fig.  1 
the  OCV— SoC  characteristics  and  the  hysteresis  between  charging 
and  discharging  is  depicted. 

SoC  estimation  is  typically  based  on  a  nonlinear  model  using 
Kalman  filter  theory.  The  nonlinear  model  describes  the  dynamic 
behaviour  of  the  terminal  voltage  in  dependance  of  the  charge/ 
discharge  current  /(f)  and  other  factors  like  e.g.  SoC.  The  SoC 
observer  is  based  on  a  combination  of  the  (relative)  SoC  model  (1) 
and  the  terminal  voltage  model,  see  Fig.  2.  Thus,  the  SoC  correction 
is  obtained  from  a  comparison  of  the  actual  terminal  voltage  U(t)  to 
the  output  of  the  model. 

In  combination  with  the  proposed  LMN  different  observer 
structures  can  be  used,  where  the  architecture  and  interpretability 
of  the  local  models  as  a  local  linearisation  of  the  process  helps  to 
reduce  the  model/observer  complexity.  These  observer  structures 
include  e.g.  the  extended  Kalman  filter,  the  interacting  multiple 
models  (I MM)  algorithm  and  the  fuzzy  observer. 

The  idea  of  the  extended  Kalman  filter  is  to  apply  conventional 
Kalman  filtering  to  a  nonlinear  system.  The  filter  gain  is  computed 
using  the  local  Jacobian  of  the  nonlinear  model.  The  main  disad¬ 
vantage  of  this  strategy  is  that  the  filter  gains  cannot  be  pre¬ 
computed  since  they  depend  on  data  (local  linearisation),  see  e.g. 
Refs.  [9,6], 

The  basic  idea  of  the  interacting  multiple  models  approach  is  to 
run  a  certain  number  of  different  linear  Kalman  filters  in  parallel, 


Fig.  2.  Schematics  of  SoC  observer  architecture. 


each  corresponding  to  a  separate  model  [27].  Mixing  of  model 
parameters  as  well  as  switching  between  the  linear  models  is 
possible.  The  determination  of  the  currently  valid  model  is  realized 
by  a  detection  scheme,  utilizing  the  validity  probability  of  the  un¬ 
derlying  models  [28]. 

In  combination  with  LMNs,  the  architecture  of  the  fuzzy 
observer  (see  e.g.  Refs.  [25,29—31])  helps  to  reduce  the  computa¬ 
tional  complexity  of  the  global  nonlinear  filter.  Each  of  the  local 
models  in  an  LMN  is  a  linear  time-invariant  dynamic  system. 
Therefore,  a  local  observer  is  designed  for  each  local  linear  state 
space  model  using  standard  Kalman  filter  theory.  The  global  filter  is 
then  derived  from  a  linear  combination  of  the  local  filters. 

3.  Battery  cell  model  identification 

In  this  section  the  nonlinear  system  identification  of  battery 
cells  for  the  purpose  of  SoC  estimation  is  presented.  The  nonlinear 
model  describes  the  dynamic  behaviour  of  the  terminal  voltage  t/(t) 
in  dependance  of  the  charge/discharge  current/]/)  and  other  factors 
like  e.g.  temperature  and  SoC.  In  the  following,  the  requirements  for 
the  cell  model  for  the  application  in  HEV  development  are  formu¬ 
lated  (c.f.  Ref.  [6]): 

•  based  only  on  readily  available  signals  only  (cell  terminal 
voltage,  cell  current,  cell  external  temperature), 

•  valid  up  to  very  high  C-rates  (±20C), 

•  highly  accurate  even  under  strong  dynamic  excitation, 

•  includes  special  effects  such  as  hysteresis  and  relaxation, 

•  suitable  for  real-time  applications. 

As  an  important  prerequisite  for  nonlinear  system  identifica¬ 
tion,  a  model  based  experiment  design  procedure  for  battery  cells  is 
presented  in  Section  3.1.  In  Section  3.2,  the  LMN  architecture  is 
briefly  described  and  adapted  for  battery  cell  modelling  in  order  to 
take  into  account  hysteresis  and  relaxation.  The  proposed  concepts 
are  then  validated  using  real  measurement  data  from  a  high  power 
Lithium  Ion  cell. 

3.1.  Optimal  experiment  design 

An  important  prerequisite  for  data  based  modelling  approaches 
is  a  suitable  experiment  design.  In  general,  the  target  of  the  DoE  is 
to  generate  informative  data  while  the  experimentation  effort  is 
reduced.  Thus,  the  experiment  design  can  be  understood  as  a 
compromise  between  experimentation  effort,  reliability  and  accu¬ 
racy  under  the  specific  loads  and  environmental  conditions. 

In  the  context  of  cell  modelling,  the  whole  operating  range  (cell 
current,  SoC)  of  the  cell  has  to  be  covered  and  the  model  has  to  be 
accurate  for  a  highly  dynamic  excitation  (high  C-rates).  One  of  the 
main  challenges  regarding  the  experiment  design,  is  that  the  SoC 
excitation  (and  the  operating  range)  directly  depends  on  the  exci¬ 
tation  signal  of  the  cell  current,  c.f.  Equation  (1).  In  nonlinear  system 
identification,  amplitude  modulated  pseudorandom  binary  signals 
(APRBS)  [19]  are  widely  used  in  order  to  track  the  nonlinearities  of 
the  underlying  process.  This  strategy  is  no  longer  feasible  for  the 
excitation  of  batteiy  cells  since  the  dependence  between  cell  current 
and  SoC  is  not  taken  into  account  which  leads  to  insufficient  coverage 
of  the  operating  range  of  the  battery.  Consequently,  the  experiment 
design  has  to  assure  that  the  variation  of  the  cell  current  on  the  one 
hand  yields  a  proper  excitation  of  the  cell  and  on  the  other  hand  that 
the  desired  SoC  operating  range  is  covered.  Additionally,  the  time  for 
battery  testing  has  to  be  taken  into  account  in  the  optimisation 
procedure.  In  this  work,  a  model  based  experiment  design  approach 
is  presented  which  is  especially  suited  for  cell  model  identification. 
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3. 1.1.  Model  based  DoE 

In  general,  model  based  design  of  experiments  is  a  methodology 
for  the  generation  of  excitation  signals,  where  the  signal  is  adjusted 
to  a  specific  prior  process  model  in  a  way  that  the  model  parame¬ 
ters  can  be  estimated  from  measured  data  with  minimal  variance, 
see  e.g.  Refs.  [17,32].  The  optimisation  can  be  based  on  linear  (e.g. 
obtained  from  preliminary  tests)  or  nonlinear  dynamic  models 
(from  a  similar  battery  cell).  In  the  present  work,  a  linear  dynamic 
cell  model,  which  was  identified  from  preliminary  tests  (step  re¬ 
sponses),  was  used  for  optimal  experiment  design.  However,  the 
results  presented  later  in  this  section  indicate  that  even  such  a 
simple  prior  process  model  helps  to  improve  the  model  quality 
significantly. 

In  order  to  maximise  the  information  of  the  excitation  signal, 
the  Fisher  information  matrix  (FIM)  is  used  [33].  In  statistical  terms 
the  Fisher  information  matrix  X  is  a  tool  to  measure  the  information 
content  of  the  excitation  signal.  The  FIM  is  a  real  valued  and  sym¬ 
metric  matrix,  whose  inverse  yields  a  lower  bound  for  the 
parameter  covariance  matrix.  The  definition  of  the  FIM  is  based  on 
the  partial  derivative  of  the  model  output  with  respect  to  the  model 
parameters,  c.f.  Ref.  [32]: 


1  N 

*  =  iZ 


9y(k ,  6)9y(k,  6y 


96 


96 


(2) 


For  the  optimisation  of  the  input  signal,  a  scalar  criterion  J  is 
required.  In  this  context,  the  D-optimality  criterion  is  widely  used 
where  the  determinant  of  the  FIM  is  subject  to  maximization: 


/fim  =  det(X). 


(3) 


The  optimisation  of  the  cell  current  excitation  signal  is  based  on 
the  following  predefined  conditions: 

•  SoC  operating  range:  To  assure  a  proper  coverage  of  the  SoC 
operating  range  (and  to  observe  SoC  limits),  predefined  SoC 
levels  have  to  be  included  in  the  experiment  design. 

•  Cell  current  excitation :  To  obtain  a  proper  distribution  of  the  cell 
current  excitation  predefined  current  levels  are  applied.  Thus, 
even  when  the  optimisation  is  based  on  a  linear  cell  model,  a 
good  coverage  of  the  input  space  is  obtained. 

•  Hysteresis  and  relaxation:  In  order  to  capture  the  voltage  hys¬ 
teresis  and  relaxation  behaviour  of  the  battery  cell,  the  charge/ 
discharge  pulses  are  followed  by  zero  current  phases. 

Thus,  the  excitation  signal  is  obtained  from  a  certain  number  of 
design  points,  each  characterised  by  the  absolute  value  of  their  cell 
current  and  SoC  level.  The  target  of  the  optimisation  procedure  is  to 
determine  an  optimal  sequential  arrangement  of  these  points.  For 
each  design  point,  the  sign  of  the  cell  current  (charging  or  dis¬ 
charging)  and  the  duration  of  the  charge/discharge  pulse  is  ob¬ 
tained  from  the  associated/desired  SoC  level  and  the  previous  SoC. 
Hence,  the  testing  time  directly  depends  on  the  sequential 
arrangement  of  the  SoC  levels  and  the  associated  current  levels.  As 
a  consequence,  the  duration  of  the  testbed  run  is  also  subject  to 
optimisation  (in  addition  to  the  information  content  of  the  signal). 

The  optimisation  is  focused  on  the  appropriate  sequential 
arrangement  of  the  SoC  levels  and  a  proper  assignment  of  the 
current  levels  using  the  simulated  annealing  method.  Similar  to  Ref. 
[34]  the  search  begins  with  an  initial  design  and  proceeds  through 
examination  of  a  sequence  of  designs,  each  generated  as  a  pertur¬ 
bation  of  the  preceding  one.  Thus,  a  new  candidate  state  is  obtained 
from  a  random  exchange  of  the  current  and  SoC  levels,  respectively. 
The  probability  of  making  a  transition  from  the  current  state  to  the 
new  state  is  specified  by  an  acceptance  probability  function. 


In  order  to  optimise  the  information  content  of  the  signal  while 
the  testing  time  is  kept  to  a  minimum,  the  following  performance 
function  is  minimised 


i  ,/FIM.init  ,  ,  ?opt 

/opt  =  a- - +  (1  -  a) 


/fim 


(4) 


opt 


where  /piM.init  denotes  the  determinant  of  the  Fisher  information 
matrix  of  the  initial  design  and  Tjnit  and  ropt  are  the  duration  of 
testbed  run  of  the  initial  and  the  optimised  design.  The  choice  of 
the  design  parameter  0  <  a  <  lin  Equation  (4)  is  based  on  the 
tradeoff  between  accuracy  and  measurement  effort:  When  a  is 
increased  the  information  content  is  increased,  while  for  a  =  0  only 
the  testing  time  is  reduced.  Note  that,  on  the  one  hand,  an 
increasing  number  of  data  also  leads  to  an  improvement  of  Equa¬ 
tion  (3).  On  the  other  hand,  the  number  of  data  cannot  be  arbitrarily 
decreased  since  the  predefined  SoC  and  current  levels  have  to  be 
included  in  the  design. 

3.1.2.  Experimental  results 

The  proposed  concepts  presented  in  this  paper  are  validated  by 
means  of  a  Lithium  Ion  cell.  Therefore,  two  training  data  records 
(obtained  from  the  initial  experiment  design  and  the  optimised 
input  signal)  and  a  validation  dataset  were  recorded.  In  this  section, 
measured  input  and  output  data  and  the  results  from  the  optimi¬ 
sation  in  Section  3.1.1  obtained  from  an  initial  design  and  the 
associated  optimised  design  are  presented. 

The  complete  initial  data  record  is  shown  in  Fig.  3,  where  the 
initial  value  of  “Charge  [Ah]”  corresponds  to  80%  SoC.  The  nominal 
capacity  of  the  power  cell  is  2  Ah  and  the  nominal  voltage  3.3  V, 
respectively.  Accordingly,  the  SoC  operating  range  is  between  20% 
and  80%.  Note  that  the  choice  of  the  initial  sequential  arrangement 
of  the  SoC  levels  plays  an  important  role  for  the  testing  time  Tjnit  of 
the  initial  dataset. 

The  initial  design  was  then  used  as  a  starting  point  for  the 
optimisation  of  Equation  (4)  using  the  simulated  annealing  algo¬ 
rithm.  Here,  the  initial  and  final  SoC  were  fixed  and  thus  not  subject 
to  optimisation.  In  Fig.  4  the  improvement  of  the  determinant  of  the 
Fisher  information  matrix  and  the  reduction  of  the  measurement 
time  are  depicted.  Obviously,  the  determinant  Equation  (3)  is 


4000  6000  8000 

time  [s] 

Fig.  3.  Training  data  record:  initial  DoE. 
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Fig.  4.  Optimisation  of  the  excitation  signal  using  simulated  annealing. 

increased  by  a  factor  of  five  while  the  testing  time  is  reduced  by 
about  7%. 

The  measured  data  record  obtained  using  the  optimised  input 
signal  is  presented  in  Fig.  5.  Again  the  SoC  operating  range  is  be¬ 
tween  20%  and  80%  and  the  predefined  current  levels  are  realised. 

3.2.  Nonlinear  system  identification 


strategy  allows  to  capture  the  highly  nonlinear  dynamic  complexity 
in  a  computationally  efficient  way. 

Each  local  model  of  the  LMN  —  indicated  by  subscript  i  —  con¬ 
sists  of  two  parts:  The  validity  function  <Pj(x(k))  and  its  model 
parameter  vector  ft,.  Thereby,  #,■  defines  the  region  of  validity  of  the 
i-th  local  model. 

The  local  estimate  for  the  output  is  obtained  by 

yt{k)  =  xT(k)Ou  (5) 

where  xT(k)  denotes  the  regressor  vector.  In  dynamic  system 
identification,  the  regressor  vector  x(k)  comprises  past  system  in¬ 
puts  and  outputs. 

All  local  estimations  yj(fc)  are  used  to  form  the  global  model 
output  y(k)  by  weighted  aggregation 

M 

y(k)  =  2^$i(k)yi(k),  (6) 

i  =  1 

where 


In  this  section,  the  identification  of  the  Lithium  Ion  power  cell 
from  the  data  presented  in  Section  3.1  is  presented.  The  input  into 
the  model  are  the  cell  current  I  and  its  past  values  up  to  an  order  of 
four,  the  state  of  charge  SoC  (k  -  1 )  and  the  model  output  is  the  cell 
terminal  voltage  U(k).  For  the  feasibility  study  presented  in  this 
paper,  the  cell  temperature  is  kept  constant. 

The  architecture  and  training  algorithm  of  the  proposed  LMN  is 
shortly  reviewed,  see  also  Ref.  [23].  Additionally,  the  adaptation 
and  extension  of  LMN  structures  for  cell  identification  (in  order  to 
take  into  account  hysteresis  and  relaxation)  is  addressed.  Simula¬ 
tion  results  using  measurement  data  from  the  Lithium  Ion  cell 
demonstrate  the  performance  of  the  proposed  model  architecture. 

3.2.1.  Local  model  network  architecture 

The  LMN  interpolates  between  different  local  models,  each  valid 
in  a  certain  region  of  the  input  space.  Thus  the  battery  cell  model  is 
based  on  a  partitioning  into  several  local  operating  regimes,  rep¬ 
resented  by  the  dominant  influence  e.g.  SoC,  temperature,  etc.  This 


0i(k)  =  $,(*(10)  (7) 

and  M  denotes  the  number  of  local  linear  models.  Thereby,  the 
elements  in  x(k)  span  the  so-called  partition  space  and  are  chosen 
on  the  basis  of  prior  knowledge  about  the  process  and  the  expected 
structure  of  its  nonlinearities.  Thus,  the  dimension  of  the  partition 
space  and  furthermore  the  complexity  of  the  optimisation  problem 
can  thus  be  reduced  dramatically,  see  also  Ref.  [19], 

The  computation  of  the  validity  functions  <P,-  is  based  on  a  lo¬ 
gistic  discriminant  tree.  In  Fig.  6  a  model  tree  with  three  local 
models  is  depicted.  Each  node  corresponds  to  a  split  of  the  partition 
space  into  two  parts  and  the  free  ends  of  the  branches  represent  the 
actual  local  models  with  their  parameter  vector  6 ,•  and  their  validity 
functions  4>i.  The  overall  nonlinear  model  thus  comprises  JVf  local 
models  and  M  -  1  nodes  which  determine  their  regions  of  validity. 

For  the  representation  of  the  discriminant  function  in  the  d-th 
node  a  logistic  sigmoid  activation  function  is  chosen,  c.f.  Ref.  [35]: 

'  1+exp< -«,(*,>)))  «8> 

with 

ad(x(k))  =  [l  xr(k)]^dd°J.  (9) 

~  T 

Here,  \j/d  =  [t/'dt  V'dp  ]  denotes  the  weight  vector  and  \j/do  is 
called  bias  term.  The  discriminant  functions  ipd  are  used  to  calculate 
the  validity  functions  <?>,•,  c.f.  Ref.  [36], 

The  validity  functions  for  the  layout  in  Fig.  6  are  obtained  by 


Fig.  5.  Training  data  record:  optimal  model  based  DoE. 


Fig.  6.  Logistic  discriminant  tree. 
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=  <Pt<P2, 

(10) 

=  1  -  <P\, 

(11) 

=  <Pl(l  -  <P2)- 

(12) 

An  incremental  model  construction  allows  to  gradually  increase 
the  complexity  of  the  local  model  network:  When  the  number  of 
local  models  M  is  increased  by  one,  the  worst  local  model  (indexed 
by  /)  of  the  logistic  discriminant  tree  in  Fig.  6  is  replaced  by  a  new 
node  and  two  adjoining  local  models  are  appended,  see  Fig.  7.  On 
the  one  hand  this  strategy  allows  a  proper  initialisation  of  the  new 
model  parameters  while  on  the  other  hand  the  computational 
demand  is  low. 

In  order  to  reduce  the  computational  demand,  only  the  weight 
vector  of  the  new  node  ^m-i  is  optimised  while  all  other  weight 
vectors  are  retained. 

An  important  topic  in  nonlinear  system  identification  with  local 
model  networks  is  the  interpretability  of  the  local  models,  see  e.g. 
Refs.  [37—39].  Only  consequent  parameters  obtained  by  local  esti¬ 
mation  (weighted  least  squares)  allow  the  interpretation  of  the 
consequent  parameters  as  a  local  linearisation  of  the  nonlinear 
system  [39].  In  order  to  obtain  an  interpretable  model  structure  for 
the  proposed  LMN,  the  following  cost  function  is  chosen  for  the 
optimisation  of  the  two  new  local  models  and  the  weight  vector  of 
the  new  node: 

1  N 

J  =  ^i(fc)l>(fc)  -yi(fc)]2  +  <Mfc)[y(k)  -yM(k)]2 

k  =  1 

=  l  Yet(k'>  +  eM(,<'>-  H3) 

k  =  1 


1.  The  values  for  the  resistance  and  capacitance  of  a  cell  model  are 
functions  of  the  direction  of  current  and  SoC,  [7[.  Thus,  separate 
models  for  charge  and  discharge  are  enforced  to  approximate 
the  voltage  hysteresis.  The  architecture  of  LMNs  allows  an  easy 
integration  of  this  dependency  using  a  user-defined  pre-par¬ 
titioning  based  on  the  current  direction  sign(i). 

2.  In  Refs.  [6],  large  time  constants  which  describe  the  relaxation 
effect  are  implemented  as  a  low-pass  filter  on  /(/<).  Accordingly, 
an  additional  input 

v  =  filt(/)  (16) 

for  the  LMN  is  used  in  this  work.  The  filt-operator  in  Equation  (16) 
represents  a  low-pass  filter  with  zero  dc  gain.  Here,  zero  dc  gain 
ensures  that  the  influence  of  the  additional  input  vanishes  after  a 
rest  period  so  that  the  model  output  converges  to  OCV,  c.f.  Ref.  [6[. 

3.3.  Modelling  results 

The  performance  of  the  LMN  is  highlighted  using  real  measure¬ 
ment  data  from  the  Lithium  Ion  power  cell.  The  training  of  the  LMN 
using  the  algorithm  presented  in  Section  3.2.1  results  in  an  LMN 
with  20  local  linear  models.  In  order  to  prove  the  generalisation 
capabilities  of  the  proposed  LMN  validation  data  were  recorded,  see 
Fig.  8.  Here,  the  current  profile  was  taken  from  Ref.  [7]  and  adjusted 
in  a  way  that  the  desired  SoC  operating  range  is  covered.  Again,  the 
initial  value  of  “Charge  [Ah]”  corresponds  to  80%  SoC. 

In  Fig.  9  a  comparison  of  measured  and  simulated  model  output 
at  validation  is  depicted.  It  is  clearly  visible  that  the  LMN  accurately 
describes  the  nonlinear  behaviour  of  the  cell  while  the  SoC  is  varied 
in  a  wide  range. 

In  Table  1  the  performance  of  the  models  obtained  from  initial 
and  the  optimised  DoE  is  indicated  by  means  of  the  R2  statistics 


where 

e,(k)  =  ]/^i(k)ly(k)  -  yj(k)} 

and 


(14) 


r2  =  1  Zk=i(y(i<)-yQ<))2 
£k=i  (y(i<)-y)2  ’ 

where  y  is  the  mean  value  of  all  y(fc). 


(17) 


eM (fc)  =  rJ<PM{k)\y{k)  -yM(k)]  (15) 

are  the  weighted  local  prediction  errors  of  the  adjoining  models  l 
and  M,  respectively.  In  Equation  (13)  N  defines  the  number  of 
training  data. 

Note  that  the  particular  choice  of  Equation  (13)  aims  at  a  locally 
weighted  least  squares  optimisation  of  the  consequent  parameters 
so  that  local  interpretability  is  conserved. 

3.2.2.  Extended  model  structure 

LMNs  offer  a  versatile  structure  for  the  identification  of  nonlinear 
dynamic  systems.  Due  to  the  transparency  of  the  LMN  architecture, 
the  incorporation  of  prior  knowledge  is  easily  possible. 

In  order  to  take  into  account  special  effects  such  as  hysteresis 
and  relaxation,  the  structure  of  the  LMN  is  adapted/extended: 


Fig.  7.  Incremental  tree  construction. 


Fig.  8.  Validation  data  record. 
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Fig.  9.  Comparison  of  simulated  model  output  and  measured  terminal  voltage 
(validation). 
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4.  Fuzzy  observer  based  SoC  estimation 

The  architecture  and  interpretability  of  LMNs  (see  Section 
3.2.1)  provides  a  framework  for  automated  nonlinear  observer 
design.  Especially,  the  use  of  a  fuzzy  observer  is  beneficial  in 
combination  with  LMNs:  Each  of  the  local  models  in  an  LMN  is  a 
linear  time-invariant  dynamic  system.  Therefore,  a  local  observer 
is  designed  for  each  local  linear  state  space  model  using  standard 
Kalman  filter  theory.  The  global  filter  is  then  derived  from  a 
linear  combination  of  the  local  filters  [25].  Thus,  the  local  ob¬ 
servers  are  time-invariant  which  greatly  reduces  the  computa¬ 
tional  complexity  of  the  global  filter  compared  to  the  widely  used 
extended  Kalman  filter.  Another  important  advantage  of  the  fuzzy 
observer  architecture  is  that  the  stability  analysis  of  the  nonlinear 
observer  is  possible  based  on  Lyapunov  stability  theory,  see  e.g. 
Refs.  [40,41], 

In  this  section,  the  fuzzy  observer  is  shortly  reviewed  and 
simulation  results  from  the  Lithium  Ion  cell  are  presented. 

4.1.  Fuzzy  observer  design 


and  the  input  vector  which  contains  current  and  past  elements  of 
the  inputs 


u(k) 


Kk) 

i(k  - 1) 
v(k) 

v{k-  1) 

v(k  -  m) 

1 


(21) 


the  state  equation  of  the  ith  local  model  is  obtained 

M 

z(k)  =  Y,  ^i(k-\){Aiz(k-\)+Biu(k)}.  (22) 

i=i 

The  system  matrix  Equation  (19)  and  the  input  matrix  Equation 
(20)  comprise  the  elements  of  the  numerator  and  denominator 
polynomials  and  the  affine  term  C;  of  the  local  affine  models. 

With  C  =  [10.  ..00]  the  global  model  output  is  obtained 


For  nonlinear  dynamic  systems  described  by  local  model  net¬ 
works,  a  fuzzy  observer  can  be  designed  to  estimate  the  state 
vector.  The  nonlinear  observer  design  involves  representation  of 
the  nonlinear  system  as  local  linear  state  space  model  [30].  Using 
standard  Kalman  filter  theory,  steady-state  filters  can  be  designed 
for  each  local  model  of  the  LMN  and  the  global  filter  is  obtained 
from  a  linear  combination  of  the  local  filters. 

The  SoC  observer  is  based  on  a  combination  of  the  (relative)  SoC 
model  Equation  (1 )  and  the  terminal  voltage  model  (see  Section  3.3). 
Thus,  an  augmented  state  space  formulation  of  the  nonlinear  model 
is  required.  Using  the  relation  SOC(k)  =  SOC(k  -  1)  +  Ts/Cni(/<)  the 
augmented  state  vector  Equation  (18)  also  comprises  SoC(k  -  1) 
which  was  originally  used  as  a  model  input  for  the  training  of  the 
LMN.  The  state  vector  is  defined  as 


Table  1 

Model  performance  at  validation  data. 


Initial  DoE 

Optimised  DoE 

R2 

0.98452 

0.99271 

y(k)  =  Cz(k)  (23) 

Using  these  relations,  a  local  steady-state  Kalman  filter  with  the 
gain  matrix  /(,  can  be  designed  for  each  local  model.  The  overall 
state  estimation  is  then  defined  as  a  weighted  sum  of  the  individual 
linear  observers.  The  state  estimate  z(k)  is  given  by  the  following 
equations: 

zj(k)  =  A(z(k  -  1 )  +  BjU(k)  (24) 

and 

M 

m  =  Y  0i(k  -  V&(k)  +  i<i\yO<)  -  y(m  (25) 

i=l 

where 

M 

y(k)  =  Y  *i(k  -l)Cz?(k)  (26) 

2=1 

The  Kalman  gain  matrix  Ki  is  obtained  from 
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I<i  =  AjPj  CT  (cpJ  CT  +  Rt)  1  (27) 

where  P,  is  the  solution  of  the  discrete-time  algebraic  Riccati 
equation  (DARE) 

AjPtAj  -  Pi  -  AtPiC T  (cPiCT  +  R)  CPtAj  +  Q  =  0  (28) 

Here,  Q  and  R  denote  the  covariance  matrices  of  the  measure¬ 
ment  noise  v  and  of  the  process  noise  p,  respectively.  In  Kalman 
filter  theory,  the  covariance  of  both  the  process  and  measurement 
noise  are  assumed  to  be  known.  However,  the  covariance  matrices 
are  often  used  as  a  design  parameter  since  in  practice  the  covari¬ 
ance  values  are  not  precisely  known.  In  the  present  application,  the 
covariance  matrix  of  the  measurement  noise  R,  can  be  estimated 
from  the  knowledge  of  the  battery  terminal  voltage.  Similar  to  Ref. 
[9],  Qis  obtained  from  experimental  studies  under  the  simplifying 
assumption  that  there  is  no  correlation  between  the  elements  of  /i. 
In  the  next  section,  SoC  estimation  results  presented  where  the 
covariance  value  of  SoC  is  used  as  a  tuning  parameter  for  the 
nonlinear  SoC  observer. 

4.2.  SoC  estimation  results 

The  performance  of  the  fuzzy  observer  for  SoC  estimation  is 
demonstrated  by  means  of  the  validation  data  record  (c.f.  Section 
3.3).  Now,  the  SoC  is  assumed  to  be  unknown  and  the  initial  state  of 
SoC  is  chosen  at  random. 

The  upper  plot  in  Fig.  10  shows  a  comparison  of  true  and  esti¬ 
mated  SoC  of  the  Lithium  Ion  cell.  The  estimated  SoC  converges  to 
the  correct  value  even  though  the  initialisation  is  chosen  outside  of 
the  operating  area  of  the  training  data.  Obviously,  the  LMN  based 
fuzzy  observer  accurately  estimates  the  unknown  SoC  and  an 
estimation  error  smaller  than  ±3%  is  obtained. 

In  Fig.  11  the  weighting  of  SoC  is  increased  in  order  to  obtain  a 
faster  convergence  after  a  poor  initialisation.  Thus,  on  the  one  hand 
a  fast  SoC  error  correction  is  obtained  while  on  the  other  hand,  the 
width  of  the  SoC  estimation  error  band  due  to  modelling  errors  and 
measurement  noise  is  slightly  increased. 


- meas. 

-  est. 


- meas. 

—  est. 


Fig.  11.  Generalisation  data:  SoC  estimation/correction  after  random  initialization. 

5.  Conclusion 

The  work  described  in  the  present  paper  shows  an  integrated 
methodology  for  experiment  design,  nonlinear  system  identifica¬ 
tion  and  nonlinear  observer  design  of  a  Lithium  Ion  power  cell  for 
the  purpose  of  SoC  estimation.  Model  based  design  of  experiments 
for  battery  cells  is  proposed.  The  results  indicate  that  even  a  simple 
prior  process  model  helps  to  improve  the  model  quality  significantly 
while  the  testing  time  is  reduced  and  the  coverage  of  the  SoC 
operating  range  is  ensured.  The  local  model  network  construction 
and  architecture  is  shortly  reviewed  and  adapted  for  battery 
modelling  in  order  to  take  into  account  hysteresis  and  relaxation 
effects.  Results  obtained  from  real  measurement  data  show  the 
excellent  generalisation  capabilities  of  this  approach.  Additionally, 
the  architecture  and  interpretability  of  the  LMN  provides  a  frame¬ 
work  for  fuzzy  observer  design.  A  local  observer  is  designed  for  each 
local  linear  state  space  model  using  standard  Kalman  filter  theory 
and  the  global  filter  is  then  derived  from  a  linear  combination  of  the 
local  filters.  Thus,  the  computational  complexity  of  the  global  filter  is 
greatly  reduced  compared  to  the  widely  used  extended  Kalman 
filter.  Again,  the  performance  of  the  proposed  fuzzy  observer  is 
demonstrated  by  means  of  real  measurement  data. 
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